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Abstract 

The leading mass dependence of the wave function phase is calculated in the 
presence of gravitational interactions. The conditions under which this phase 
contains terms depending on both the square of the mass and the gravitational 
constant are determined. The observability of such terms is briefly discussed. 



a. Introduction There has been a series of publications discussing the mass dependence 
of the gravitational effects on the phase of the wave function and, due to a divergence 

in the results, a certain amount of controversy has resulted. Some results 0] indicate the 
presence of a gravitationally-induced phase proportional to the square of the masses, while 
other calculations show that such terms are in fact absent 0], and that the leading depen- 
dence in the gravitational contributions to the phase is proportional to the fourth power of 
the mass. In this short note I will attempt to resolve this (apparent) contradiction. I will 
first consider the case of neutrinos propagating in an arbitrary gravitational field and then 
generalize to other particles 



2 



b. Neutrinos in a gravitational field Using a WKB approximation it is possible to derive 
the effective Hamiltonian for neutrinos propagating in a non-trivial metric. Following the 
procedure described in detail in the relevant expression is [| 



^abcd _ 2^cib^cdef PePf 



^2 



75 + lrn\ (1) 



where Latin indices refer to the local Lorentz frame, Greek indices to the global coordinate 
frame, and denote the tetrads, p is the momentum of a null geodesic which the neutrino 
wave packets follow, within the WKB approximation, when m = 0. I will denote by A the 
affine parameter along this geodesic 

The vector p denotes the time component of p with respect to a comoving reference frame. [| 
The Hamiltonian H^s generates translations in A. 

Though the pure gravitational terms in (|l|) can generate interesting effects, here I will 
concentrate on the contributions generated by the term iv?/2. It is clear from the above 
expression that each mass eigenstate acquires a phase equal to 

= -^m^A. (3) 

There are subleading contributions to i^cfr which depend on the mass @], but these are 
of order m/R where R denotes the distance scale of the metric, and can be ignored for 
m > 10-"eV. 

c. Other particles The same expression for the mass dependence of the phase can be 
obtained using the following argument. Within the WKB approximation the wave function 
of a general quantum system takes the form 



^The Hamiltonian has units of mass^ due to the choice of evolution parameter, see below 
^More precisely: denoting by v^, A = 1,2,3 three independent (local) solutions to the geodesic 
deviation equation ||5|, p is the component of p orthogonal to the u'^. 
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^ = e^^x, (4) 

where S denotes the classical action and x ^ slowly-varying amplitude. When the particle 
moves in a non-trivial metric background S satisfies the following Hamilton- Jacobi equa- 
tion 

g'^'d^Sd^S = m\ (5) 

(I ignored all other interactions and considered a single mass eigenstate). The action is 
generally covariant so that the corresponding phase is unambiguously determined (within 
the semiclassical approximation) 

For small m the action takes the form 

S = So + m'Si + • • • , (6) 

which, when substituted into @ yields 

g'^'p^Pu = 0, p^d.S, = -i (7) 

where 

= -g^^d^So = — , (8) 

is the tangent vector to the null geodesic. It then follows that p^d^Si = dSi/dX, so that 
Si = — which gives the same contribution as (|^). It follows that this result is general 

provided m is small compared to the particle's momentum. 

d. Observability of the gravitationally induced phase I will consider a wave packet that 
can be decomposed into a coherent superposition of mass eigenstates. Using standard argu- 
ments the above expressions imply that, within the semiclassical approximation, the phase 
5*0 + may lead to oscillations between components of different mass. This will occur 

provided (a) the amplitude of the mass components of the state have a sufficient overlap at 
the observation point and, (b) the phase difference between them is sufficiently large (^vr). 
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I will first review briefly the conditions for under which there is significant amplitude overlap 
and then discuss the effects of the phase difference. 

If the mass components are created in the same localized space-time region A then, for 
long distances and within the WKB approximation, they will overlap at the observation 
region B only if the corresponding amplitudes in momentum space are centered at different 
momenta. In this case the phase difference is simply /\m?{\B — A^)/2. 

Alternatively we can assume that the initial mass amplitudes are localized in the same 
region of momentum space. Then, in order for them to arrive simultaneously at B (so that 
there is significant amplitude overlap), they must originate at different space-time locations 
A and A', which may be space or tiime-like separated (the creation region then contains A 
and A' and is sufficiently large to allow a small momentum uncertainty). In this case the 
phase difference is Sq{A) — Sq{A') + /S.rn?{XB — Ayi)/2, where Sq{A) — Sq{A') is also of order 
Am^. This is the scenario described in 0. 

Leaving aside the issue of which of the two scenarios above is more natural, it is clear that 
the expressions for the phase difference will be different in each case because the assumed 
initial conditions are different |Q. 

Even if interference effects are present, an independent question is whether the phase 
difference between two mass eigenstates has a strong dependence on the gravitational field. 
I will consider for simplicity the first case described above, where the phase for each mass 
component is given by (^ and depends only on the geodesic length A. In order to determine 
the dependence on the gravitational force this expression should be re-written in terms of 
a set of physical observables, but such a set is no unique and different choices may exhibit 
varying dependences on the gravitational constant. 

The expression for will depend on the geodesic parameters such as the energy and 
angular momentum. If, for example, the interference effects are measured by a freely-falling 
observer, these should be written in terms of the corresponding quantities measured by the 
observer (e.g. the energy- momentum vector equals e^^p^). In addition, this phase will depend 
on the physical distance form the source to the observer. This distance can be defined as 
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half the time it takes for hght to go from A to B and back to A (in units where the speed 
of hght is one); exphcitly 

iAB = / dXJ{gij - goigQj/gQo)p'p^ = E / ——= (9) 
Ja Ja y/-goo 

(Latin indices denote the spatial components) the second expression is valid for null geodesies 

in a static space time (the integrals are evaluated along the geodesic). Note however, that 

Eab is not the only possible quantity that can be associated with the physical distance 

from the source to the observer. Other possibilities are in general available and, as will be 

illustrated below, the form of are strongly dependent on the choice made. 

There is one final "experimental" issue concerning the observability of the phase differ- 
ence. In many situations the distance to the source is only approximately known and the 
location of creation region is known only with limited accuracy. In such cases the gravita- 
tional effects will be distinguishable only if the phase difference due to this uncertainty can 
be neglected or otherwise taken into account. 

e. Gravitational dependence of the quantum phase in a Kerr metric. The expression 
(§) contains a dependence on the gravitational interactions due to the non-trivial expression 
of A in terms of the coordinates. As an illustrative example I will consider the case of Kerr 
metric, 

ds-'=(l-!lA de - ^dr^ - p^d9^ -(r^ + a^ + ^ sin^ o] sin^ Od^'' + ^ sin^ Od^dt, 
\ J A \ P^ J P^ 

(10) 

where 

A = - rrg + p'^ = r"^ + cos^ Q. (11) 

The gravitational radius equals 2MG where M is the mass of the black hole and G 
Newton's constant. The black hole angular momentum equals aM. 

In this case the Hamilton- Jacobi equation (^) with m = is separable, 

5o = -Et + Wr{r) + We{d) + L0, (12) 
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where E denotes the energy and L the (azimuthal) angular momentum for the particle (the 
z axis is along the the black hole rotation axis). Both E and L are constant, the addition 
constant of the motion K is related to the total angular momentum. Substituting, 



\ dd j \ sin 6'/ V dr 

Using (H) it follows that 



A=/ , ^^^^ (14) 

'[{r^ + a^)E - aLf - KA] 



with p and A defined in (|TT]) where 9 (buried in A) is assumed to be expressed in terms of 
r using 

which also follows from (^. The flat-space limit of ( p!4D corresponds to = a = and 
K = L^, in this case 



Aflat space = ^Vr2-(L/E)2, (16) 

which is a standard result. Note that A has units of mass^^, so that H^s has units of mass^ 
as noted earlier. 

To lowest order in r^, a, K — I find, 

A^l7^^T(Z7^_(^-^^) + 2^^--^^-.+.... (17) 

2E^^r^ - {L/Ef 

In the limit of zero angular momentum 0fC = L^ = O,a = O this expression retains a term 
proportional to r^, but this term is actually an irrelevant constant, rg/{2E). It follows that 
for weak gravitational fields and small masses (^ becomes 



$M = _!!lj r2_iL/Ey + -^'[iK-L-)+2ELa] _ m\{r - ^r^ - {LjEf) ^ 
2^^ 4^3 Jr2 - {LjEf AEJr^ - {L/Ef 



(18) 



a is the black-hole angular momentum per unit mass 
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where I explicitly subtracted the above-mentioned constant. It follows that when 
written in terms ofr, does receive gravitational contributions, but these are proportional to 
an angular momentum (squared). This dependence on the angular momentum significantly 
suppresses the magnitude of these gravitational terms. In a region very distant from the 
black hole, r ^ the above expression reduces to (see also 0) 

^ = - W |l + ^ + ■ ■ j ' » L/E). (19) 

For zero orbital angular momentum there is still a contribution to (^) proportional to r^, 
but it is also proportional to the square of the black-hole angular momentum. For example, 
for a geodesic in the 6 = 7t/2 plane 

$(-) = -|^(^l + ^ + ---); {L = 0e = n/2). (20) 

The expressions (0,^) contain no explicit gravitational contribution when the geodesies 
and the black hole have zero angular momentum. However when expressed in terms of the 
invariant length (|^) the phase (p!8D becomes 



- ^arcsinh(^£/L) + ■ • 



(21) 



2E 

which does contain a term ~ r^m^ In £ which persists even in the L limit (note that 
in this limit there is a term proportional to InL, but it is constant and will cancel when 
considering phase differences). Yet, despite appearances, (pTSD and (^) are equivalent. 

The calculation presented in was done for radial geodesies in a Schwarzschild metric 
and corresponds to the result ( |T^) when K = L = a = 0, in which case there is no explicit 
gravitational terms present in the phase. The same expression corresponds to the results 
in provided the quantity J dL/r used in these references (where L represents a distance 
and should not be confused with the angular momentum as used in this paper) is identified 
with 2 In i. It is worth emphasizing that both i and r have reasonable physical interpretations 
as measures of the distance to the source ^ and both reduce tO the usual flat-space distance 



^For example, in the case a = 0, r can be defined in terms of the circumference of circles around 
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in the limit of zero gravitational interactions. 

In order to extract the gravitational term from these expressions one must subtract the 
non-gravitational contribution. This, however, is an ambiguous quantity. For example, 
an observer defining the fiat space contribution as = —m^r/{2E) will find a small 

gravitational contribution of order (ang. momentum) ^/(ri?)^ to (^. A second observer 
might alternatively define ^^^^l = —vn^ij {2E) and find a much more significant gravitational 
contribution ~ r^lnf to Equivalently, an experiment done in fiat space over a distance 
d will generate a phase —m'^d/{2E), in comparing this to a similar measurement in curved 
space, one must determine what quantity (r, i, etc.) is to be used in lieu of d, and the 
expression for depends significantly on the choice made. In this sense the results 

presented in and are, in fact, consistent. 

As mentioned in paragraph |^ the the observability of these effects is a separate issue, 
which I discuss briefly in paragraph |^ below. Here I merely note that the phase (in terms 
of r) is of order {rn^ / E){r + r|/r) when L ~ ^gE:, the second term is absent for the case 
of zero angular momentum. In terms of ^ and the phase contains a term ~ {m?rg/ E) ln£, 
even when the angular momenta vanish. It is possible, however, that these simple results 
are consequences of the high degree of symmetry of the metric (|TD|). This possibility is 
investigated below. 

/. General spherically symmetric metric The metric for this configuration is using 
polar coordinates, 

ds^ = e^de - e^dr^ - r\d9^ + sin^ 9d^^), (22) 

where u and 7 depend on r and t only. The geodesies of this metric lie on a plane which I 
take as the 6 = 7r/2; then (|^) with m = is again separable and its solution takes the form 
So = -Et + L(j) + Wr{r) where = e^iE^e-" - L^r^) so that 



l/^^-^ p3) 



the origin while the interpretation of £ was provided in paragraph |^. 
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which can in principle be solved for t(A), r(A) and inverted to obtain A = A(r). For the 
interesting case of time-independent metric this gives 



$M = / dr ^ (24) 

It appears possible for this to retain a gravitational contribution even when L = (radial 
geodesies). To determine whether this is in fact the case I will consider the case of an ideal 
fluid of energy density p and pressure p both of which vanish for r > R. The solution to 
Einstein's equations gives 

pr 

7 = ln[l — 2m(r)G'/r], m{r) = / ATTu'^p{u)du 

Jo 

V = -2G / du— 7n-^7^^ (25) 

Jr u[u — 2m{u)G\ 

where G denotes Newton's constant; for r > R u and 7 reduce to their Schwarzschild 
expressions. To first order in G I find 



roc 

1 + AnG / du u[p{u) + p{u)] (26) 

Jr 



which, when substituted into (|24D with L = gives the phase for radial motion. If the 
initial point of the geodesic rj satisfies ri < R and final point lies beyond R then a simple 
estimate gives "l*^™^ ~ [m? /2E){r + cvg) where denotes the gravitational radius of the 
matter distribution, and c is a numerical constant ^O(l) that depends on the detailed form 



of p and p. Note that, just as in (|1^,^), the leading gravitational contribution to the phase 
depends on the initial point and the details of the metric, but not on r; if expressed in terms 
of £, however, the phase again acquires an m? In C. contribution. 

g. Weak gravitational interactions For the case of a general weak metric it is possible 
to obtain a rather simple expression for the phase <l>(™). Expanding Sq (cf. (P)) in powers 
of the metric perturbation, 

g,u = ri>.. + K.; 5o = 4°^ + 4'^ + ---; S'^^=0{h), (27) 

(where rj denotes the flat-space metric) and defining 
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(0) 



(28) 



I find — and 



(29) 



where indices are raised and lowered using the flat metric 77. 
It proves convenient to choose coordinates such that 

p^ ^ E{1,1,0,0), 

then the equation for Sq^^ becomes 



(1) 



^^ = -2(^00 + 2/101 + ^11), 



(30) 



(31) 



so that 



So ^ -p ■ X - J dx+{hoo + 2/ioi + hu) + 



(32) 



Prom this and the Hamilton- Jacobi equation I find 

dx^^ 



dX 

and, in particular, 



^ pM + ^Edi" J dx+{hoo + 2/ioi + /ill) - h'^'^P. + 



c/x''^ 1 f 

^E- -Edi / cix+(/ioo + 2/ioi + /in) + E{hu + /iio) + • • • ■ 
aA 2 J 

For a time-independent metric the above expression simplifies to 

dr^ 1 

^^E--E{hoo-hn), 



whence XE = x + ^ J dx{hoo — hu) and 



2£; 



a; + ^ J dx{hQo - hu) + 



(33) 



(34) 



(35) 



(36) 



where the phase difference must be obtained by evaluating the quantity in brackets at two 
different points on a geodesic with the line integral being along the same geodesic. 
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The gravitational terms vanish if the geodesic has zero angular momentum for the case 
of a distant (localized) matter distribution since in this case |^ /iqo = —"^g/^ and hij = 



dependence on the angular momenta will be quadratic. This result is not accidental: linear 
terms are absent since the phase is a scalar and the angular momentum is a pseudovector. 
In terms of the invariant length 



which again has a linear+logarithmic dependence on £. 

h. Discussion The leading mass dependence of the phase of a general wave function 
is, within the WKB approximation, proportional to the affine parameter along the geodesic 
followed by the wave packets and may contain explicit gravitationally-induced contributions. 
The expression for the gravitationally-induced phase depends on the initial conditions as- 
sumed for the system and on the choice of physical variables. This last point is particularly 
delicate since the separation of into gravitational and fiat-space contributions is sensi- 
tive to the quantity used as physical distance. The results presented in are consistent 
with each other provided these points are taken into account. 

In all interesting cases, however, the uncertainty in the distances to regions of strong 
gravitational field are too large for these effects to be unambiguously isolated experimentally. 
In terms of the radial coordinate r (defined at large distances using the Schwarzschild metric) 
and for the case of non-zero angular momentum, ~ —{m'^/2E){r + r^jri) where 



denotes the initial point of the geodesic and tq is a characteristic length of the problem 
(such as the geodesic impact parameter or the gravitational angular momentum per unit 
mass) such that the angular momenta are of order r^E. In all cases considered tq < < rj 
so that ~ —{w? l2E){r + r^). These gravitational effects (for the case where the flat- 
space phase is defined as —m?r/2E) can be isolated experimentally only when r is known 
to an accuracy better than r^, and this is not usually attainable. 

In terms of the invariant length (. the phase takes the form ~ {rn?/2E)[l + 




Cartesian coordinates). Moreover, since "l'^'"-* is independent of /loi, the 




(37) 
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{rg/2)\n.{l/rg)\ with a more significant dependence on Vg. Nonetheless the precision in I 
required to extract the gravitational contribution (where the flat-space phase is now deflned 
as w?l/2E) is still not available. For example for a lOM© supernova 10^ light years away, 
I should be known better to a precision better than 2, 000 km or 10~^^ %; for a lO^M© 
supermassive quasar 10^ light years away, £ should be known to a precision better than 
10^^ km or 10~^° %. Despite the logarithmic enhancement the current experimental sen- 
sitivity precludes an accurate discrimination of the gravitational contributions for realistic 
situations. 
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